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ABSTRACT 


This thesis evaluates the so-called Liouville-Green 
approximaticn to the solution of variable-coefficient 
Lanchester-type equations for combat between two homogeneous 
forces, When compared to the form of the exact solutions, 
this approximation is in terms of "elementary" functions. 
Two specific forms of attrition-rate coefficients are 
considered, allowing for different maximum effective ranges 
of the two ofposing weapon systems. These coefficients might 
be used to model a constant-speed attack against a static 
defensive position. Le 1s shown that for these 
attrition-rate coefficients, the Liouville-Green 
approximation is not consistently reliable for predicting 
force levels, and yields exact results only under certain 
restrictive conditions. Furthermore it was found that 
methodology is not presently available to accurately predict 
from Liouville's normal form the error which will be 
incurred by invoking the approximation in a specific 


Situation. 
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I. INTRODUCTION 


Although it is well recognized that combat is a complex 
random process, a model that provides "closed-form" 
solutions can often be an aid to the military analyst by 
explictly portraying interrelationships between combat 
parameters. Deterministic differential equations are 
frequently used to model the combat process, and vary 
greatly in their detail and level of complexity. The 
Military analyst must be able to recognize significant 
variables and their relationships. An idealized model is 
useful for giving a "first cut" portrayal of such 
relationships. In particular, one can often obtain an 
analytical sclution to such models. An analytical solution 
that can be expressed in terms of "elementary" functions can 
often provide insights into the combat process and identify 


factors affecting battle outcome. 


The use of differential equations to model warfare was 
pioneered by F. W. Lanchester in 1914 [Ref 8]. This original 
work was quite idealized and assumed that the fire 
effectivenesses were constant throughout the battle. Many 
extensions to Lanchester's classic formulation have 
Subsequently been developed, to include such factors as 
mobility of forces and fire effectiveness that varies during 
the battle. Taylor and Brown {Ref 13] have developed a 
mathematical theory for solving Lanchester-type equations 
With these temporal variations in fire effectiveness. The 
fire effectiveness of forces is represented by what are now 
called Lanchester attrition-rate coefficients. Taylor and 
Brown give results for two specific forms of the 
coefficients: (1) effectiveness of each side's fire 
proportional to a power of time, and (2) effectiveness of 
each side's fire linear with time but a nonconstant ratio of 


attrition- rate coefficients. 
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Taylor and Brown present accurate numerical solutions 
to equations with the above attrition-rate coefficients. 
These solutions are, in general, stated in terms of infinite 
series. The solution method, however, is rather complicated 
and involved the development of new mathematical functions 
alien to the military analyst. Exact results are certainly 
beneficial, but the analyst also needs a solution in terms 
of "familiar" functions to permit parametric analysis and 
‘simplify recognition of significant variables. To provide 
such a solution form, Taylor [Ref 14] has suggested the use 
of the so~called Liouville-Green approximation, which can be 
expressed in terms of relatively simple functions. For an 
approximation to be useful in predicting force levels, 
however, it must be reasonably "close" to the exact 


solution. 


This paper provides a numerical evaluation of the 
Liouville-green approximation for the two forms of 
attrition-rate coefficients mentioned above. A comparison is 
made between the approximation in each of these cases and 
the numerical results given by Taylor and Brown. Numerical 
examples and analytic considerations will be utilized to 
investigate the accuracy of the approximation in these 


cases. 


The organization of this thesis is as follows. We first 
review Lanchester's classic model of combat between two 
homogeneous forces. These results are then extended to the 
variable attrition-rate coefficient case, and the specific 
forms to be studied are discussed. The representation of 
the solution to these variable-coefficient equations in 
terms of so-called general Lanchester functions is 
presented. Then Liouville's normal form and approximation 
are introduced, evaluated, and discussed for the two 


specific forms of attrition-rate coefficients. 





II. LANCHESTER'S CLASSIC FORMULATION 


F. W. Lanchester was an English aeronautical engineer 
who first developed a mathematical formulation for combat 
between two conflicting forces. He published this theory in 
a series of articles in the British journal Engineering in 
1914 [Ref 8]. Lanchester's objectives were to provide 
insights into "modern" combat and to quantitatively justify 
the military principle of concentration. He stressed that 
advances in technology had brought long-range delivery 
capability of modern weapons which allowed for concentration 
of firepower. Lanchester hypothesized that under these 
modern conditions, combat between two military forces could 
be modelled by: 


dx/dt ~ay with x (t=0) 


li 
ra 
< 


is (1) 


dy/dt 
Y/ 0 


- bx With y (t=0) 
where x(t) and y(t) are the numbers of X and Y forces at 
time t, and t=0 denotes the time at which the battle begins. 
These equations are valid only for x and y greater than 
zero. The parameters a and b are non-negative constants that 
have come to be referred to as Lanchester attrition-rate 
coefficients. These coefficients represent each side's fire 
effectiveness or firepower. The most common set of 
circumstances under which these equations have been 
hypothesized to apply 1s that both sides use aimed fire and 
target acquisition times are constant [Ref 16]. 


From (1) Lanchester deduced his classic square law 
2 2 Z 2 
Oe at Cie) e) S aa - y (t))- (2), 


which has the important implication that a force can 


Significantly reduce its casualties by initiaily committing 
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more forces to the battle. The square law also shows that in 
this simple model the Y force wins (ie., X is annihilated) 


if and only if 
1 a < Va/b. (3) 


From (1) Lanchester also derived the now well-known 
results for the time history of the force levels 


alt) = x coshVab t - y va/b Saowaivabe tr 
vt ils coshVab t - x Vbsa SinhVab t , (4) 


The victory-prediction condition (3) was given for a 
"fight to the finish". As H. K. Weiss [Ref 15] points out, 
most battles terminate before complete annihilation of one 
force occurs. Different models of battle termination have 
been proposed, among them that combat ends when either of 
two given "breakpoint" force ratios is reached. Taylor has 
Shown that (3) is necessary and sufficient for Y to "win" a 


fixed-force-ratio-breakpoint battle [Ref 14]. 


Lanchester's original formulation (1) involved many 
implicit assumptions that limit its usefulness for modelling 
real combat. The forces were assumed to be homogeneous, 
there were no replacements or withdrawals, and no movement 
Of forces was considered. Setting the attrition-rate 
coefficients a and b in (1) as constants implies that the 
fire effectiveness per firer of each side is constant 
throughout the battle. 
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III. VARIABLE ATTRITION-RATE COEFFICIENTS 


When modelling mobile weapon systems and movement of 
forces, one recognizes that the assumption of constant fire 
effectiveness is open to question [Ref 2]. T£ fire 
effectiveness is permitted to vary during the battle, 


Lanchester's formulation (1) takes the forn: 


dx/dt a(t) y with x (t=0) 


il 
a) 
&% 


i 
ns 
ba) 
-— 

O1 
— 


dy/fdt = -b(t) x with y (t=0) 


where a(t) and b(t) are time-dependent attrition-rate 
coefficients. The form of the attrition-rate coefficients 
depends on such variables as force separation(range of 
battle), tactical posture of the targets, firing rate and 


rate of target acquisition [Ref 5]. 


New operations research techniques for predicting these 
attrition-rate coefficients were developed in the 1960's. A 
Significant contribution in this area was the development of 
a methodology by S. Bonder [Ref 2 and 4] for prediction of 
coefficients based on weapon system performance data. 
Another important development was the method of G. Clark for 
the maximum likelihood estimation of these coefficients 
based on Monte Carlo simulation output data [Ref 6]. The 
work of these and others has facilitated the use of nmodels 


such as (5) in current combat studies. 


To effectively model today's highly mobile warfare, the 
form of the attrition-rate coefficients in (5) must 
accomodate current tactics and battle conditions. 
Specifically, the coefficients Must permit the fire 
effectiveness of one or both sides to be non-zero at the 
Start of battle {to model, for example, an ambush). The 
coefficients must also allow for the opposing weapon systems 


to have different maximum effective ranges. Since mobile 
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forces close with each other in combat, their fire 
effectiveness should be permitted to increase or decrease as 
a function of range or time. The attrition-rate coefficients 
examined in this paper have the flexibility to nodel these 
various Situations. Motivated by previous work by Bonder 
and Farrell [Ref 5], Taylor and Brown have considered the 


following attrition-rate coefficients; 
m n 
a(t) = k (t+#C) and b(t) = Kk (tecea) j (6) 
a 


The specific cases of (6) to be considered here will be that 
in which A=0, c20, 


mi n 
a(t) = k_(t#C) and b(t) = k (t#C) : (7) 


and that in which m=n=1, A>0O, C20, 


a(t) = k (t+C) and b(t) = Soe (8) 
a 


The parameter A will be referred to asthe "offset 
parametertt since it permits modelling of combat between 
forces whose weapon systems have different maximum effective 
ranges. The parameter A also allows the ratio of 
attrition-rate coefficients to be non-constant. Koopman [see 
Ref 10] first observed that if the ratio of attrition-rate 
coefficients is constant, the time solution to (5) is no 
more complicated than the constant-coefficient solution. 


That is, when 


a(t) = k h(t) and b(t) = kK b(t) « (9) 
a 


where h(t) is the common time-dependent factor in both 


attrition-rate coefficients, the solution is 


x(t) = x cosh W(t) ~ y Vk /kK sinh W(t), 
0 O' a b 
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where W(t) = Vk re Frcs) ds. Koopman also found that the 
a 0 
constant ratio of coefficients again yields a square-law 


relationship(see also Ref 16). 


The parameter C in (7) and (8) allows for the initial 
fire effectiveness of each side to be non-zero. This permits 
modelling of battles which begin within the naxinun 
effective ranges of the opposing weapon systems. The 
parameters m and n allow the fire effectiveness of each side 
to increase or decrease with time. With m,n<0, the fire 
effectivenesses decrease. This situation might model a 
battle in which combatants take cover and improve their 
positions as the battle progresses. With m,n>0, the 
effectivenesses increase, such as might be the case in which 


two forces are closing with each other. 


To model an engagement in which forces close with one 
another, it is beneficial to relate their instantaneous fire 
effectiveness to the range separating them. As this range 
decreases, one might expect fire effectiveness to increase. 
The attrition-rate coefficients (6) can be expressed in 
terms of range between forces and the maxinum effective 


ranges of the opposing weapon systems: 


dx/fdt 


m 
SOE} ay e= -42 (1-7 ae yy, 

0 a 
dy/fdt = -b(r) x = -b (I-r/R) x, (10) 


These equations have been used by S. Bonder [Refs 1 and 3] 
to model a constant-speed attack on a static defensive 
position. FPigures 1 and 2 show sample relationships between 
the attrition-rate coefficients and range between opposing 
forces for various input parameters. In these equations, fr 


is the range between the forces, and R and ae are the 
a 
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maximum effective ranges of the Y and X weapon systems, 
respectively. The constants ar, and , represent the fire 
effectiveness of each side at r=0. Range and time are 
related by ESR Vite where Be denotes the opening range of 


battle, and v is the constant attack speed. With this 
relationship, the attrition-rate coefficients in (10) can be 
related to (7) and (8). For (7) we have: 


m n 
k =a (vV/R k = B (V/R C = (R -R V 
a Pe a 5 io Can e ie ne “(11) 
where R = R 2R, while for (8) we have: 
a b 0) 
k =a v/R k = BV/R A = (R -R )/v, C= (R -B )/¥ 
a 0 a ; b 0 B’ b oa ’ a 0O : 
where R > R 2R. (12) 
b a 0 
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IV. A STANDARD FORM FOR THE VARIABLE-COEFFICIENT SOLUTION 


The variable attrition-rate coefficient equations (5) 
yield the force level equations 
2 2 
d x/fdt - [1/fa(t) da/dt] dx/sdt - a(t) b(t) x = 0, (13) 
With initial conditions 


x(t=0) = x and 1fa(t) dx/dt =-y , 
and Similarly 


2 2 
d y/fdt - [1/b(t) db/dt] dyydt - a(t)b(t) y = 0, (14) 
with initial conditions 


t=0) = d 1/b(t) dy/dt = -xX . 
ae la BOSC SGD ion a x 


The solution to (13) is given by 


Cy = 70 SK oe CGE, PS 
x (t) ere tee At) (5D) 


where pee el sha | denotes a fundamental system of 
solutions to (13), and S and ee, are constants determined by 


the initial conditions [{see, for example, Ref 7]. The 


Y-force level has similar form. The functions SOE Se) G 
Se and ee are referred to by Taylor and Brown as 
General Lanchester Functions. Their properties are 
Summarized in Table I. The functions ee! and es are 
Similar to the hyperbolic cosine, while a) and co) are 
Similar to the hyperbolic sine. The General Lanchester 


Functions may be constructed by successive approximations or 


by infinite series methods. 
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A. APPLICATICN TO POWER ATTRITION-RATE COEFFICIENTS 


When this solution method is applied to the case of 
power attrition-rate coefficients (7), the solution (15) may 


be written as 


matt) =x [0 (C) u (t+C) - V (C)v (t+C) ] - 
0 n n m,n 


’ Ry, ¢ ’ 
k /k u C)v tC) tal v C 16 
YJ a7 b [ en ) Bast ) ae oth ) ], (16) 
where u EMS , JU and V acre referred to by Taylor as 
M,N m,n m,n m,n 


Power Lanchester Functions. They have the properties shown 
in Table Ii. For computational convenience, utilization is 
made of the Lanchester-Clifford-Schlafli (LCS) functions 


. 2k K-] 
F (x) = 5(x/2)” / (k 1G+v) 1, 
Vv k=0 j=0 
2k+1 k 
G (X) = 2 (x/2) 4 (k Il (jtv) J. 
k=0 j=0 


The Power Lanchester Functions may be expressed in terms of 


the LCS functions through the following relationships: 


u (t) = F (S(t)) 
n,n gq 

(m-n) /2 
v cp eee G (S (t)) 
m,n p 
U (t) = F (S(t)) 
n,o p 

~ 2 

VS gee ae Ge 
m,n g 


where S(t) =\/ (k ty (kt) t / ((mtnt2)/2). Using the tLcs 
a 


functions, (16) may now be written as 
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(m—-n) /2 
ceeey = Xe {od (0)) F (J (e)) = (1+t/C) G (J(0))G (J(t)) 
0 p q q p 


(m—n) /2 
- Y va (t=0) /b (t=0) [ (1+t/C) We (J baybes se cc)) = 


G (J(0))F (J (t)) }. 17 
5 : 7) 
B. APPLICATION TO OFFSET LINEAR COEFFICIENTS 


Similar techniques are utilized to apply this general 
theory to the case of the offset linear attrition-rate 
coefficients (8). With these coefficients, a fundamental 


system of solutions to (13) is given by 


oe (t) =f (t+C), “ (t) =g (ttC), Y, (t) =F (t+C), Y, (t) =G (t+tC) 


where 
fore) 2n k k 4n-k 
f(t) = EE {VK k /2) /s(2ny!} TB At 
n=0 a b k=QON 
ee) 2n+¢1 nk k Unte2~-k 
g(t) = r{Wk k, 7/2) /(2atiy!} FC At 
n=0 a D k=O fl 
co 2n n k k 4n-k 
Pit), = 2 10K 72) fuego) i} sp A Ee 
n=0 aD k=0 8 
00 2n+1 n+lk k 4n+2-k 
G(t) = © ¢Wk k 7/2) 7Zizortt) yo Bea. t : 
n=O ab k=0 8 
is k k k k 
The coefficients B,C, D, and & Satisfy complicated 
n n n n 


recurrence relations [see Taylor and Brown (Ref 13) j}. The 


functions £(t), g(t), F(t), and G(t) are referred to as 
Offset Linear Lanchester Functions. As in the case of tne 
power attrition-rate coefficients, use is made of the 
following Auxiliary Offset Linear Lanchester Functions for 


computaticnal convenience: 
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00 2n nk k 
h(p,u) = %{p /(2n)!} 2B p 

n=0 k=Q2 

oO 2n+1 n k k 
w(Pp,u) = 2& {p /(2n+1):1} FTC u 

n=0 k=0 7 

o 2n nk k 
H(p,u) = L£{p /(2n)'!} 2D u 

n=0 k=Q A 

co «= 2+ 1 n+l k 
Wtp,u) = 2 {p /(2nt+i!} EBu. 

n=0 k=0 2 


Utilizing these Auxillary Lanchester Functions, the X~-force 


level solution (15) may be written as 
E(t) = x CHT (0) 65 (0)) (T(t) oJ (EI) — W(T(0) +5 (0)) 

w(J (t) 5 (t)) J 
- y VK /&, Ch (3 (0) +5 (0)) w(F(t) ,5(t)) - w(5(0) ,4(0)) 


h(J (t),J3(t)) ], (18) 
2 
where j(t) = A/(t+C) and J(t) =Vk aM (EFC) Jas 
a 


The equations (17) and (18) are solutions to the 
X-force level equation (13) for the case of power and offset 
linear attrition-rate coefficients, respectively. This 
sclution method of Taylor and Brown represents a 
considerable extension of Lanchester theory in the area of 
Variable attrition-rate coefficients. Previous analytic 
results had been available only in the more limited case of 
the opening range of battle equal to the minimum cf the 
maximum effective ranges of the two weapon systems. That is, 
the initial fire effectiveness of at least one side was 


restricted to be zero [see for example Ref 12]. 
While the time solutions (17) and (18) =appear 


formidable, their implementation is Simplified if a digital 


computer is available. As noted earlier in this paper, 
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however, these solutions are complex and involve new 
mathematical functions. This makes it difficult for the 
military analyst to perceive significant relationships among 
battle variables, and renders parametric analysis difficult 
and time-consuming. This provides motivation for the 


approximation technique to follow. 
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V. REDUCTION TO LIOUVILLE'S NORMAL FORM AND APPROXIMATION 


ee pee a 





The argument ab (referred to by, among others, Taylor 
and Brown [Ref 13] as the "intensity of combat") in the 
constant coefficient result (4) provides motivation for the 


transformation [Ref 14] 
t 
z= sfVa(s)b(s) ds, (19) 
ft 
where z(t=0) is denoted Zs Applying this transformation to 


(13) and (14) yields 


2 2 
d x/dz + {1/2 d/dz ln{b(t)/a(t) ]}dx/dz - x = 0, (20) 


With initial conditions 


1/2 
X(Z=Z ) = xX and {{ b(t) /a(t) ] f dx7az} = -y , 
0 0 Z=Z 0 


and 


2 2 
d y/fdz + {1/2 d/dz Infa(t)/b(t) J} dyydz - y = 0, (23) 


with initial conditions 


We 
y(z=z _) = y and {f[a(t)/b(t) ] dy/dz} =i Xue 
8) 0 Z=Z ¢ 


Taylor and Brown have shown that the X force level 
equation (13) may be transformed into a linear second order 
differential equation with constant coefficients if and only 
1f (20) is a constant coefficient equation. This is the case 
1£f and only if 


(1/Va(t)b(t)) d/dt In{a(t)/b(t) ] = CONSTANT. 


Equation (20) shows the significance of the parameters 


a(t)/b(t) (the relative effectiveness of the two weapon 
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systems), and the intensity OL combat Va (t) b(t) in 
determining battle outcome. This 1S a generalization of the 
well-known constant coefficient results to the case of 
temporal variations in fire effectiveness. The significance 
of these parameters may be _ seen nore explictly by 


transformation to Liouville's normal form (see Ref 14). 
The transformation 
1/4 
X(Z) = SCN Se / ee 72 ey | ‘ ae) 


may be applied to (20) and yields the normal form with the 


first derivative of the dependent variable removed 


2 2 
d X/faz - {1 + F(z) ]X = 0, (ee) 


with initial conditions 


Ola = ad axysd = =-y b - 
(Z omy an an JaAZ{Z aa ue Thy ; oa 
where 
rt -1/4 
F(z) = P (z) / P(z) P(z) = [R(t) ] (24) 
R(t) = a(t) /b{t) e(t) = 1/4 Va(t) b(t) d/dt LnR(t). 


(25) 
2 2 
In (24), the notation P''(z) denotes d P/dz . Equation (23) 


is Liouville's normal form [see page 23 of Ref 9}. 


For the Y force level, the transformation 


1/4 
y(z) = pA CUE 4) PNG ED LN aC 


yields the normal forn 


Pa Z 
d Y/fdz - [1+ G(z) ]Y = 0, (26) 


With initial conditions 
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Y (Z=Z = and dY/dz(z=z ~x Vb fa + iy 
( 9? Yo {AZ ( o! , 0” 0 Yo 


0 


where 


1/4 
G(z) = Q**(z) / Q(z) and Q(z) = ([R(t)J = 1/P(z). 


Expressing the transformed force-level equation (23) in 


the simple forn 


2 2 
d X/dz - X = Q(z)X, 


variation of parameters may be used to obtain the solution 
to (23) as 


x = h (Z- = V b + Xx inh (z~ 
{Z) my Nie ae CY, ie ; hae Sinh (z aa) 
+ {F(4) sinh (z-4) X (5) aj. (27) 


Zz 
This equation may be expressed in terms of the original time 


variable t and dependent variable x as 


1/4 C 
x(t) = Fale od we) Wn | Gat icra fois) 5 (eras) 
i 


“Ly Va /b, + x .@, Jsimh ( /va (s) b (s) ds) } 
0 


1/4 1/74 
+ Ca(t)/b(t) J d/ds{1/(a(s)9(s))d/ds(b(s)/a(s)] 3 
sinh( f¥Va(sj)b(j) dj) X(s) ds. (28) 


Ss 
If the relative effectiveness varies slowly, then the 


relationship (24) would indicate that F(z)<<1 so that one 
could drop the integral term in (27) yielding 


X(Z) = a cosa (ZeZ)! - Cy Va sb, + SAT ae ee oe 


which is the Liouville-Green approximation [Ref 11]. The 
approximation may be expressed in terms of the original 


variables t and x as 
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1/74 C 
x(t) = See See ‘ ee a (s) b (s) ds) 
t 0) 
= DSN + x25 Jsinh ( /¥a (s) b (s) ds) } (30) 
0) 
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VI. PPLICATION TO POWER AND OFFSET LINEAR COEFFICIENTS 


The Liouville-Green approximation may now be applied to 
the constant-speed attack nodel With attrition-rate 
coefficients (7). For the reader's convenience, these 


coefficients are 
Nn n 
a(t) = k (t+tC) and b(t) = k (t+C) P (7) 
a 


This model assumes that the maximum effective ranges of the 
two opposing weapon systems are equal (no offset). The 
opening range of battie is taken to be within this maxinun 
effective range (i.e., C>0). Applying the Liouville-Green 
approximation (29), the time solution to the X¥ force level 
equation in terms of the original time variable t and 


dependent variable X becomes 


= 4 £ 
x(t) = Meth ayes (x cosh /Va (s)b(s) ds 
0 


= Ze -(mtnt+2) /2 
- {y Vk /k °@ i +(x (m-n)/4 k k )C aaa ] 
QO a ob 0 ab 


c 
Seiya Ye CS) CIS) say) 


where 


t 
f¥a(s)b(s)ds = (Vk k / (ment 2) /2) 
0 a 


Ge at oh aero taal aa 


For future reference, equation (24) becomes in this 
case 


Cee 
F(z) = ({ (mn) (3mtnt4) ] / (4(mtnte) Zz ], 32) 


or in terms of the original time variable t and range 


parameters (11): 
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ret} = [RR (m-n) (3m+n+4) ] / 
a b 
16 b aL, =O t+(R -R L gat A 
a ve - ve a 
[ 0 9 | ) (t+ ( ; ale ) J ( 388) 


where vel = constant velocity of attack. 


We can also apply the Liouville-Green approximation to 


the case of the offset Linear attrition-rate coefficients 


a(t) = k (t+C) and b(t) = Sy (eefeuls) (8) 
a 


In this case, the time solution to the X force level 


equation beccmes 


1/4 t 
x(t) = [ (1+ (A/C) / (1+ (AZ (t#C))) J t (x cosh /Ya(s)b(s) ds 
0 


-{(y Vk /k )/ (1+ (A/C)) + (x (A/SC))/ 
| 0 a »b 0) 
2 372 ees 
(4 k kC (1+ (A/C) ) ) JSinhSYa(s)b(s) ds}, (34) 
a 0 
where 
: 2 2 
{Va (5) b (s} ds = (Aw k k 78) { (1+ (2 (t+C) 7a) ) (14+ (2(t+C) /A)) -1 
a 


2 
—- [nf 1+ (2(t#C)/A) + (14+(2(tt+C) /A)) -1)} 


For future reference, equation (24) becomes in this 


case 


F(t) = R RB (R -R ){12(tt+(B -R )/vel) + 7((R -R )/vel)] / 
aB Boa a QO B a 


3 S| 3 
(16a b vel (t+(R -B )/vel) (t+ (BR -R ) /vel) Dim (Sue 
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VII. DISCUSSION 


A computer program was developed in FORTRAN to produce 
plots of the xX force level versus range between opposing 
forces for the constant-speed attack model described 
earlier. Both the power and offset linear attrition-rate 
coefficient cases (equations (31) and (34)) were considered. 
These approximate solutions were compared to the exact 
results given by Taylor and Brown {Ref 13] to evaluate the 
accuracy of the Liouville-Green approximation to equations 


with these coefficients. 


A. POWER ATTRITION-RATE COEFFICIENTS 


For the case of power attrition-rate coefficients, the 
Maximum effective ranges of the two weapon systems were set 
equal to 2000 meters, with an opening range of battle of 
1250 meters. The fire effectiveness of the two sides at zero 
range separation was Si voe X casualties/min-Y¥Y unit, and 
pe ° Y casualties/min-X unit. The velocity of attack was 
constant at five miles per hour, and the initial force 
levels were X = 10 and Y = 30. These parameter values 
correspond to those used by Tayler and Brown in their paper 
[Ref 13]. The exponents m and n in (7) were varied 
throughout the range explored in Ref 13. The results are 
Shown in Figures 3 through 10, where each figure contains 
the approximate time solution and the corresponding exact 
analytical solution obtained by Taylor and Brown. AS even a 
CUursory examination reveals, the approximation appears to 
conform to the established results only in certain cases. 
Some possible explanations for this observed deviation will 


now be considered. 


Recall that the integral term in equation (27) was 


Za 








dropped (assumed to be approximately zero) in making the 
Liouville-Green approximation. The validity of this 
assumption appears to be linked strongly to the form of 


P(z), which for the present case is 


F(t) = [RR (m-n) (3ntne4) ] / 
a 


m+n m+n+2 
[16a b (vel) ic J- (3) 


When the exponents m and n are equal, F(t)=0, and by 
considering equation (27), one can see that the 
approximation is, in fact, the exact solution. This is 
verified trivially by Pigures 3 and 4, for which the 
exponents are equal. It should be noted that in this 
Situation, the model reduces to the case of a constant ratio 


of attrition-rate coefficients. 


When the exponent m is greater than n, F(t) will always 
be greater than zero, and from (27) it follows that the 
exact soluticn is greater than the approximation. This 
Situation is illustrated by Figures 6 and 10. The converse 
also holds. When m is less than n, F(t) is always less than 
zero, and the approximation 1S greater than the exact 


solution. This is shown in Figures 5 and 9. 


Again examining (33), the parameter|ja b |, which is the 


product of the opposing forces! fire effectivenesses at r=0, 


appears to be significant. As ja | increases, F(t) 
decreases, making the integral term in (27) smaller, so one 


would expect the approximation to be "closer" to the exact 
solution. This did occur in the cases considered (see for 


example Pigure 7). 
B. OF PSET LINEAR COEFFICIENTS 
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For the case of the offset linear attrition-rate 
coefficients (8), the maximum effective range of the Y force 


weapon system (R) was set equal to 1500 meters, with the 
a 
rest of the parameter values equal to those in the power 


coefficient case. The maximum effective range of the X force 


weapon system (RB) waS varied throughout the range 
considered in Ref 13. The results are Shown in Figures 11 


through 15 where, as before, each figure compares the 
corresponding solutions. AS was the case with the power 
coefficients, the solutions are seen to conform only in 


certain cases. 
For the linear coefficients, it was shown that F(z) 
could be written in terms of the time variable t as 
P(t) = RB R (R -R )EI2Z(t+(R -R )/fvel) + 7((R -R )/vel) j / 
aB Boa a QO Boa 
3 : 3 
{16a b vel (t+(R -R )/vel) (tt(R -R )/vel) }. (35) 
00 a 0 b 0 
When the saximum effective ranges of the opposing weapon 
systems (R and a) are equal, F(t) = 0, and the integral 
a 
term in (27) is zero. This again reduces to the constant 


attrition-rate-ratio case as can be seen from (8) with A=0. 
Thus one would again expect agreement of the exact and 
approximate solutions. This was found to be true as shown 
by Figure 11, where both maximum effective ranges equal 1500 


meters. 


AS was the case for the power coefficients, the 


parameter Jato appears in the denominator of F(t), hence 


the two solutions should be closer as ab | increases. This 


Z9 





was found to be the case, as illustrated by Figure 15. 


C. ANALYTICAL CONSIDERATIONS 


The examinations of both the power and offset linear 
attrition-rate coefficient cases have centered on the 
parameter F(z) in the solution (27). While F(z) does appear 
to be the driving factor in the integral term in (27) (thus 
in the validity of the approximation), it Should be made 
clear that the other factors in the integral term can not be 
ignored. While a mathematical analysis(for example of the 
hyperbolic sine term) would be formidible, examples can be 
found which illustrate that F(z) can not always be 
considered alone. One such example is shown in Figure 8. The 
fire effectiveness of the two forces at r=0 were 


interchanged (i.e., ae and Dees’ instead of vice versa) 
in the power attrition-rate coefficient case. From the 


appropriate form of F(z), equation (33), one can see that 
this change has no effect on the numerical value of F(t) 
throughout the time range. Yet the Significant difference in 
results can be seen in Figure 8 when compared to Pigure 6. 
Currently, there is no convenient methodology to analyze the 


lmplicaticns of such discrepancies. 


Toward this end, theoretical analysis in the area of 
Olver's "Error Bounds for the Liouville-Green Approximation" 
{Ref 11] could be very beneficial. While they sound 
promising, Olver's resultS can not be directly applied to 
variable-coefficient Lanchester-type differential equations. 
If available in a convenient form, such error bounds could 
be evaluated for specific Situations to investigate 


tolerable error limits prior to invoking the approximation. 


From the above analytic considerations, one can _ see 
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that while one may be able to explain some of the error 
incurred, this is small comfort to the analyst attempting to 
make practical use of the Liouville-Green approximation. It 
has been shown that the approximation appears to yield 
credible results only in certain limited circumstances. 
Furthermore, methodology is not presently available to 
predict the error which will be incurred in a specific 
Situation, except for the trivial case in which the integral 


term in (27) goes to zefro. 
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VIII. SUMMARY 


This thesis has examined the "adequacy" of the 
so-called Liouville~Green approximation to the solution of 
variable-coefficient Lanchester-type equations of modern 
warfare. Although this approximation Only involves 
"elementary" functions and has an intuitively appealing 
forn, sige was unfortunately found that this 
Liouville-Green-Lanchester approximation is not consistently 
reliable for .estimating force levels. Specifically, the 
approximation was applied to two forms of attrition-rate 
coefficients: (I) power attrition-rate coefficients 
(modelling, for example, the same maximum effective ranges 
for the opposing weapon systems), and (II) offset linear 
attrition-rate coefficients (modelling different maxinun 
effective ranges). It was seen that the approximation was 
made by dropping the integral term in equation (27), hence 
exact results are achieved only when this term equals zero. 
This situation was found to occur only when the ratio of 
attrition-rate coefficients was constant. For the power 
attrition-rate coefficients, this implies that the exponents 
in (7) are equal, while for the offset linear coefficients 
this implies the offset parameter A in (8) 1s zero (which 
occurs when the maximum effective ranges are equal). 
Analytic considerations showed that while certain initial 
parameter combinations produced credible results from the 
approximation, there was no methodology available to predict 
the error incurred in the general case. It was found that, 
aS a general rule, the further one moves from the constant 
coefficient ratio case, the larger the error incurred by 


invoking the Liouville-Green approximation. 


a2 


_ 





~ 
- 





TABLE I. 
Properties of the General Lanchester Functions X12 Xoo Vy> Vo- 


ie dx,/dt = vk Kk, a(t)y, 


2: dx,/dt = Vig /k, a(t)y, 
S}e dy, /dt = vk /k, b(t)x, 
ie dy,/dt = Ve Tie, b(t)x, 
5. x, (ey, (t) - x (t)yj(t) =1 Ve 


- 


616 





26 


Table II. Properties of the Power Lanchester Functions 


u > Vv ee » and V : 
m,n m,n m,n m,n 
mM 
du, ,/dt ~ vig /k. [kt Ny 


Zi 
i 


™m 
/k [kt ee 


n a 
du, ,/at = Ye Te (hye vy, Ct) 

Wy n/dt = YEAMe (eye dy, g(t) 

eC Uc? = eee = 1 Wit 
eae = ew) 

Vin a (t=O) = Vian 6t=0) = 0 

Ug mot) = Ug g(t) = cosh (He ny (etl) ) 
Va git) = V(t) = sinh (fe ig cE Gat) 


rd 


Properties 6 and 7 only hold for m,n ? -l. 
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